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We  present  a  new  approach  that  incorporates  two  models  to  estimate  the  underwater 
light  field  from  remote  sensing  of  ocean  color.  The  first  employs  a  series  of  analytical, 
semi-analytical,  and  empirical  algorithms  to  retrieve  the  spectrum  of  inherent  opti¬ 
cal  properties  (IOPs),  including  the  absorption  and  the  backscattcr  coefficients,  from 
the  spectrum  of  remote  sensing  reflectance.  The  second  model  computes  the  profile  of 
photosynthetiealiy  available  radiation  ?(z)  for  a  vertically  homogeneous  water 
column  using  the  information  of  the  retrieved  IOPs  and  the  ambient  optical  environ¬ 
ment.  This  computation  is  based  on  an  improved  look-up  table  technology  that  pos¬ 
sesses  high  accuracy,  comparable  with  the  full  solution  of  the  radiative  transfer  equa¬ 
tion,  and  meets  the  computational  requirement  of  remote  sensing  application.  This 
new  approach  was  validated  by  in  situ  measurements  and  an  extensive  model-to-modcl 
comparison  with  a  wide  range  of  IOPs.  We  successfully  mapped  the  compensation 
depth  by  applying  this  new  approach  to  process  the  SeaWiFS  imagery.  This  research 
suggests  that  Eq  PAR(z)  can  he  obtained  routinely  from  ocean-color  data  and  may  have 
significant  implications  for  the  estimation  of  global  heat  and  carbon  budget. 


Keywords: 

*  Photosynthetiealiy 
available  radiation, 

■  ocean  color, 

■  radiative  transfer, 

■  ocean  optics, 

■  compensation 
depth. 


1.  Introduction 

Observing  ocean  color  from  space  is  an  important 
international  Loo!  in  resource  management  and  scientific 
investigations.  Currently,  there  are  ten  space-borne  ocean- 
color  sensors  orbiting  the  Earth  operated  by  various  space 
agencies  in  the  United  States,  Asia  and  Europe,  Another 
six  satellites  with  ocean-color  sensors  onboard  are  sched¬ 
uled  to  be  launched  in  the  next  four  years  (http:// 
www.joceg.org/sensorsjoccg.htrnl).  A  major  goal  of  these 
missions  is  to  assess  the  role  that  the  ocean  plays  in  the 
global  carbon  cycle  and  to  examine  the  factors  that  affect 
global  climate  change  ( Hooker  etaL,  1992),  To  attain  this 
goal,  the  amount,  distribution  and  productivity  of 
phytoplankton  must  be  estimated  from  space.  Many  past 
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efforts  have  focused  on  relating  radiative  signals  to  pig¬ 
ment  concentration  (Gordon  et  aL ,  1988;  O’Reilly  et  aL, 
1998;  Carder  et  aLt  1999),  Maps  of  satellite-derived  pig¬ 
ment  concentration  (e*g.,  biomass)  have  been  widely  used 
to  estimate  global  ocean  carbon  content  and  productivity 
(Longhurst,  1995;  Platt  et  at. ,  1995).  However,  rates  of 
phytoplankton  photosynthesis  are  regulated  by  many  fac¬ 
tors,  including  light  availability.  Hence,  an  accurate  esti¬ 
mation  of  ocean  primary  production  and  carbon  flux  re¬ 
quires  a  thorough  description  of  the  underwater  light  field. 
Compared  to  advanced  methods  for  retrieving  pig¬ 
ment  concentrations  from  ocean  color  (Gordon  et  al.t 
1988;  O’Reilly  et  a/.,  1998;  Carder  et  aL,  1999),  meth¬ 
ods  for  estimating  the  underwater  light  field  from  space 
have  not  progressed  very  far.  A  common  approach  is  to 
diminish  the  surface  solar  irradiance  based  on  satellite- 
derived  estimates  of  the  diffuse  attenuation  coefficient 
for  the  downwelling  planar  irradiance  KEJ (symbols  used 
in  this  paper  are  summarized  in  appendices)  at  490  nm. 
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1  able  1.  Various  values  of  the  coefficients  used  in  Etj.  (I)  taken  from  various  past  works. 
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Based  on  a  simple  linear  regression  analysis  on  the  data 
set,  Austin  and  Petzold  (1981)  proposed  the  first  empiri¬ 
cal  algorithm  of  KS(J,49Q)  of  the  form: 


0)  -  K£(f  miler  (490)  +  A 


M*t) 


(J) 


where  K£dtWalet{490)  is  the  diffuse  attenuation  coefficient 
for  pure  water  at  wavelength  490  nm,  LW(X j)  and  LJ A^) 
are  water-leaving  radiances  at  the  respective  wavelengths 
of  X{  and  A^  After  new  profiles  were  added  to  the  data 
set,  the  values  of  the  coefficients  used  in  Eq.  (1)  had  to 
be  revised  to  maintain  a  better  modeLto-data  fit.  Table  1 
lists  various  results  of  regression  analysis  from  a  series 
of  works  performed  in  the  past. 

There  are  four  major  known  sources  of  error,  how¬ 
ever,  involved  in  using  this  KEd(4 90) -based  approach. 
First,  K£d  should  be  spectrally  dependent.  Measurements 
of  Kej  at  one  spectral  band  (e.g.  490  nm)  cannot  repre¬ 
sent  the  spectrally  dependent  absorption  of  short-  and 
long-wave  components  of  light,  even  in  the  upper  10  m 
of  the  ocean  (Simpson  and  Dickey,  1981).  Second,  KEd 
should  be  a  function  of  depth  even  in  a  homogeneous 
medium  (Zaneveld,  1989).  Taking  an  average  value  of 
KEd over  the  first  optical  depth  to  propagate  the  light  field 
throughout  the  water  column  may  generate  large  errors 
at  greater  depths.  Third,  light  absorption  by  phytoplankton 
or  other  water  constituents  generally  have  no  preferen¬ 
tial  direction  (Kirk,  1994).  Current  KEd  models  only  cal¬ 
culate  the  downwelling  planar  irradiance  Ed,  yet  should 
include  the  scalar  irradiance  E0  for  heal  and  carbon  budg¬ 
ets.  Fourth,  the  underwater  light  field  is  closely  related 
to  the  sky  radiance  (diffuse)  distribution  at  the  sea  sur¬ 
face,  which  is  very  sensitive  to  solar  position  (Liu  et  al 
2002),  A  value  of  KEd  obtained  during  a  satellite  over¬ 
pass  provides  no  information  on  the  underwater  light  field 
at  other  times  during  the  day. 

Over  the  past  two  decades,  major  advances  in 
radiative  transfer  theory  have  enabled  a  realistic  numeri¬ 
cal  simulation  of  the  underwater  light  field  (Mobley  et 
aL,  1993),  provided  that  the  inherent  optical  properties 
(lOPs)  and  the  ambient  optical  environment  are  given. 
These  radiative  transfer  models,  however,  have  not  at¬ 


tempted  to  conduct  global  simulations  at  the  spatial  reso¬ 
lution  of  ocean-color  sensors,  mainly  due  to  the  signifi¬ 
cant  computational  effort  required  to  solve  the  full 
radiative  transfer  equation. 

We  present  a  new  approach  that  incorporates  two 
models  to  estimate  the  underwater  light  field  from  remote 
sensing  of  ocean  color.  The  first  employs  a  series  of  ana¬ 
lytical,  semi-analytical  and  empirical  algorithms  to  re¬ 
trieve  the  spectrum  of  inherent  optical  properties  (IGPs), 
including  the  absorption  and  the  backscatter  coefficients, 
from  the  spectrum  of  remote  sensing  reflectance  (Lee  et 
oI.t  2002).  The  second  model  computes  the  profile  of 
photosynlhetically  available  radiation  E0PAR(z)  fora  ver¬ 
tically  homogeneous  water  column  based  on  the  infor¬ 
mation  of  the  retrieved  lOPs  and  the  ambient  optical  en¬ 
vironment  (Liu  et  al. ,  2002),  This  computation  is  based 
on  an  improved  look-up  table  (LUT)  technology  that  pos¬ 
sesses  a  high  accuracy,  comparable  with  the  full  solution 
of  the  radiative  transfer  equation  and  meets  the  computa¬ 
tional  requirement  of  remote  sensing  application. 

This  new  approach  was  first  validated  against  in  situ 
measurements.  To  further  examine  the  applicability  of  this 
new  approach  for  a  large  variety  of  water  types,  an  ex¬ 
tensive  model-to-modd  comparison  was  then  performed 
for  a  wide  range  of  IOP  combinations.  The  results  show 
that  the  KEd  approach  generates  substantial  errors  under 
different  computational  conditions,  while  our  model  re¬ 
duces  the  errors  significantly.  This  new  approach  was  also 
applied  to  the  SeaWiFS  imagery  to  generate  a  map  of  the 
compensation  depth.  The  premise  assumptions  of  the  set 
of  three-variable  bio-optical  model  and  the  vertically  ho¬ 
mogeneous  water  column  have  been  carefully  discussed. 
This  research  suggests  that  EQPAJt(z)  can  be  obtained  rou¬ 
tinely  from  ocean-color  data  and  may  have  significant 
implications  for  the  estimation  of  global  heat  and  carbon 
budget, 

2,  Methods 

2,  /  Fast t  accurate  model  of  underwater  scalar  irradi¬ 
ance 

Liu  et  al.  (2002)  employed  four  strategies  to  accel¬ 
erate  the  simulation  of  the  underwater  light  field  without 
losing  accuracy  compared  with  the  latest  version  4.2  of 
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Hydrolight  (H42).  They  demonstrated  that  the  sky  radi¬ 
ance  can  be  reallocated  to  the  plane  of  the  Sun  and  de¬ 
composed  into  ten  independent  sources  of  light  incident 
from  ten  different  zenith  angles.  E0PAR(z)  can  then  be 
obtained  by  summation  of  contributions  from  these  light 
sources.  For  each  incident  light  source  with  unit  inten¬ 
sity,  they  constructed  a  look-up  table  (LUT)  for  quick 
reference  to  a  set  of  parameters  (Z?0,  Bt,  P,  B2,  Q )  that  are 
required  by  a  five-parameter  model  (McCormick,  1995) 
for  specifying  the  vertical  profile  of  the  average  cosine 
]J(z).  Finally,  Gershun’s  equation  (Gershun,  1939)  was 
used  to  convert  ft  (z)  to  £0(z),  They  also  developed  an 
empirical  approach  to  correct  the  effects  of  a  wind-rough¬ 
ened  surface  and  the  contributions  from  an  extra  CDOM 
component.  Their  model  runs  more  than  fourteen  thou¬ 
sand  times  faster  than  the  full  Hydrolight  code,  while  lim¬ 
iting  the  percentage  error  to  2.20%  and  the  maximum  er¬ 
ror  to  less  than  4.78%.  Detailed  procedures  are  illustrated 
in  the  block  of  Step  3  of  Fig.  2. 

To  construct  a  LUT  of  fi(z)r  Liu  et  al.  (2002)  em¬ 
ployed  a  set  of  bio-optical  models  to  parameterize  lOPs, 
including  absorption,  scattering,  and  the  phase  scattering 
function.  The  absorption  coefficient  a( A)  is  obtained  us¬ 
ing  the  models  of  Morel  (1991)  and  Prieur  and 
Sathyendranath  (1981): 

a(X)  =  aw(X)  +  0.06a;  (A)C/j/065 

+F  0.06  a*c (440)C/j/°  65  exp(-0.QI4(A  -  440)},  (2) 

where  aw(X)  is  the  absorption  coefficient  of  pure  water, 
Chi  is  the  chlorophy  11-a  concentration  (mg  m-3),  F speci¬ 
fies  how  CDOM  absorption  is  set  to  be  proportional  to 
chlorophyll  absorption  at  a  reference  wavelength  440  nm, 
and  ac*( A)  is  the  normalized  chlorophyll-specific  absorp¬ 
tion  coefficient.  aw(X)  is  taken  from  Pope  and  Fry  ( 1 997). 
af*(A)  is  given  by  Morel  (1988)  and  is  equal  to  unity  at  a 
reference  wavelength  440  nm  (i,e.  ar*(440)=  1).  Although 
F  is  usually  set  at  a  fixed  value  of  0.2  for  Case  1  waters 
(Prieur  and  Sathyendranath,  1981;  Morel,  1991),  it  may 
vary  within  a  range  around  0.2  (Liu  etaL,  1999).  We  de¬ 
termine  F  from  an  estimate  of  the  CDOM  absorption  co¬ 
efficient  at  440  nm  a^(440).  The  scattering  coefficient  b(X) 
is  derived  from  the  mode!  of  {Gordon  and  Morel,  1983): 

b(X)  =  6W(A)  +  0.30^  ^pjcA/0,62,  (3) 

where  bw(X)  is  the  spectral  scattering  coefficient  of  pure 
water  as  taken  from  (Morel,  1974).  The  normalized  phase 
scattering  function  $  is  the  sum  of  the  contributions  by 
pure  water  and  particles  (Mobley,  1994) 


Fig.  L  Comparison  of  derived  o(670)  (Eq,  (I))  and  original 
values  of  the  simulated  data. 


(4) 

where  ^  is  the  scattering  angle  and  hp(X)  is  the  spectral 
scattering  coefficient  of  particles,  and  Pp(yr,  A) 

are  calculated  from  the  analytic  Fournier-Forand  phase 
function  (Fournier  and  Forand,  1994)  by  specifying  the 
backscatter  fraction  of  water  BFW  and  particles  BFp ,  re¬ 
spectively,  The  value  of  BFW  is  set  to  be  0.5  (Mobley  and 
Sundman,  200 1 ), 

This  set  of  three-variable  (Chi,  F ,  BFp)  bio-optical 
models  gives  a  flexible  parameterization  of  IOPs  repre¬ 
sentative  of  a  large  variety  of  water  types.  The  current 
LUT  is  constructed  by  changing  variables  within  certain 
ranges:  A  (400-700  nm),  Chi  (0—  1 0  mg  nr3)  and  BFp 
(0.01-0.04)  (Liu  et  al.,  2002).  Note  that  there  is  no  limi¬ 
tation  on  the  range  of  F  by  employing  an  empirical  ap¬ 
proach  to  correct  the  contributions  from  an  extra  CDOM 
component. 

In  principle,  the  model  of  Liu  et  aL  (2002)  can  be 
integrated  with  any  ocean  color  algorithms.  We  present 
here  the  integration  of  their  model  and  the  quasi-analyii- 
cal  algorithm  (QAA)  (Lee  et  al.,  2002),  with  the  goal  of 
improving  the  retrievals  of  Chi ,  F  and  BF  ,  and  hence, 
the  estimation  of  E0PAR(z)  from  the  spectrum  of  remote 
sensing  reflectance  measured  above  the  surface  Rrs( A), 

2. 2  Quas i-an a lyti cal  a Igorith m 

A  full  description  of  the  QAA  is  presented  by  Lee  et 
al.  (2002).  Briefly,  the  QAA  employs  a  series  of  analyti¬ 
cal,  semi-analytical  and  empirical  algorithms  to  convert 
Rrs( X)  to  a(X)  and  bb( A).  Multiple  H42  simulations  were 
done  to  derive  the  empirical  relationship  between  a(670) 
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Fig,  2,  Processing  steps  to  convert  Rr,(X)  to  E0(z). 


and  the  ratio  of  r„(670)  and  rfJ{443)  by  varying  the  three 
variables  within  the  ranges  C/i/  (0—10  mg-nr3),  F  (0- 1 ,0) 
and  BFp  (0,005-0,05),  where  rrs  stands  for  the  remote 
sensing  reflectance  measured  just  below  the  surface.  A 
regression  analysis  gives 


o(670)  =  0.439  +  0.305 


r„(670) 


1 0.824 


Lr„(443)J 


(5) 


The  comparison  of  the  derived  a(670)  and  their  original 
values  of  simulated  data  is  shown  in  Fig.  1 .  Even  at  band 
670  nm,  where  absorption  is  mainly  dominated  by  water, 
varying  the  volume  scattering  phase  function  (VSF) 
through  BFp  has  an  effect  on  rrs>  and  hence  the  derived 
a(670).  The  simple  regression  analysis  shown  here  (Eq. 
(5))  provides  a  good  estimation  of  a(67Q).  Recent  work 
by  Lee  and  Sandidge  (2003)  employed  a  neural  network 
approach  to  relate  a  to  the  spectral  ratio  of  rrs.  Therefore, 
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Fig.  3.  Comparison  between  retrieval  of  (a)  Chi ,  (b)  F  and 
(c)  BFp  using  QAA  and  original  values  used  in  H42 
simulations.  A  total  of  100  cases  were  compared  by  ran¬ 
domly  specifying  values  of  the  solar  zenith  angle  (0°— 45°), 
percent  cloud  cover  (0-50%),  surface  wind  speed  (0-15 
m/s),  chlorophyll  concentration  (0-9  mg  m"3),  visibility  (5- 
50  km),  the  CDQM  ratio  F (0. 1 5-1.0)  and  backscatter  frac¬ 
tion  BFp  (0*0 1 5-0.025). 


the  deviations  shown  in  Fig.  1  can  be  reduced* 

QAA  was  derived  from  a  large  combination  of  lOPs 
for  a  wide  variety  of  water  types.  It  should  be  applicable 
to  Case  t  waters  of  which  the  FOPs  are  described  by  a  set 
of  three-variable  biooptical  models  (Eqs*  (2)-{4))*  The 


Table  2.  Accuracy  of  retrievals  for  simulations  shown  in  Fig. 
3. 


Index 

Chi 

F 

Pearson  correlation  coefficient 
Average  value  of  relative  error  (%) 

0.998 

5.19 

0.998 

1.73 

0.967 

6.78 

accuracy  of  QAA  in  retrieving  (Chi  Ft  BFp)  was  first 
examined  by  use  of  Hydrolight-simulated  Rrs( A)  curves. 
Detailed  procedures  are  illustrated  in  the  block  of  Step  I 
of  Fig*  2* 

A  total  of  100  cases  were  simulated  by  randomly 
specifying  values  of  seven  parameters  (see  caption  of  Fig. 
3  for  details)*  For  each  case,  H42  was  run  to  simulate  /?rr 
at  20  wavebands  in  the  PAR  range  (default  setting  in  H42 
for  simulating  SeaWiFS).  Note  that  all  processes  of  in¬ 
elastic  scattering  are  considered  and  the  fluorescence 
efficiencies  are  ihe  default  settings  in  H42.  The  simulated 
R rj(A)  were  then  used  as  input  to  QAA  for  converting 
/?„(A)  to  a(X)  and  hb(X).  The  algebraic  solutions  of  Chi 
F  and  BFp  can  . be  obtained  from  £7(412),  a(443)  and 
^(443)  (see  the  block  of  Step  2  of  Fig.  2)*  The  retrievals 
of  Chi  F  and  BFp  and  their  original  values  are  compared 
in  Fig.  3,  and  the  average  values  of  the  relative  error  for 
retrieving  Chi,  Fand  BFp  are  listed  in  Table  2* 

23  Field  data 

In  situ  data  were  collected  during  EcoHAB  cruises 
between  March  1999  and  October  2001  on  the  West 
Florida  Shelf,  and  the  Coastal  Benlhic  Optical  Properties 
(CoBOP)  experiment  from  1998  to  2000  in  late  May  to 
early  June  around  Lee  Stocking  Island  in  the  Bahamas. 
Rrs(X)  was  determined  by  correcting  for  the  surface-re¬ 
flected  skylight  and  solar  glint  from  the  above-surface 
total  reflectance,  which  is  measured  from  a  hyper-spec¬ 
tral  handheld  radiometer  (Spectrix)  developed  by  the 
University  of  South  Florida*  The  details  of  the  method 
are  given  in  Lee  et  ai  (1998).  The  vertical  profile  of  Ej(z) 
was  measured  using  a  submersible  hyperspectral  radiom¬ 
eter  that  was  also  developed  by  the  University  of  South 
Florida*  This  instrument,  together  with  two  ac-9™  ab¬ 
sorption  meters  (WET  Labs)  and  a  Hydroseat-6™ 
backscattering  meter  (HOBI  Labs)  were  deployed  on  a 
slow-drop  package  and  allowed  to  drift  away  from  the 
boat  lo  minimize  ship  shadow. 

Surface  water  samples  were  collected  at  each  sta¬ 
tion*  The  spectrophotometric  methods  consisted  of  sum¬ 
ming  the  particulate  absorption  (ap)  as  determined  by  the 
quantitative  filler  pad  method  (Bissett  et  ai ,  1997)  and 
the  absorption  of  the  0.22  fm\  filtrate  (ag)  as  measured 
using  a  10  cm  cell  in  a  Perkin  Elmer  Lambda  18  spectro¬ 
photometer  (Spec}*  The  accuracy  of  the  measurement 
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Table  3.  Descriptions  of  the  in  situ  data  used  in  the  model-to-data  comparison  (Figs.  4  and  5). 


Originator  identifier 

LSD5260I 

LK404U 

Mirirl 

Mirir2 

Location 

LSI 

(23.8 1DN.  76.06°W) 

WFS 

(26,98°NI  83 J0PW) 

(28.79*N,  89,84°W) 

(28.68*N,  89.90° W) 

Water  type 

Low  Chl-fl 

Moderate  Chi -a 

High  Chi -a 

Moderate  Chl-o 

Moderate  CDOM 

Moderate  CDOM 

High  CDOM 

High  CDOM 

Date 

2  6- May -00 

22-Apr-0I 

7 -Apr-00 

26 -Oc  1-00 

Time  (GMT) 

13:57 

18:03 

14:10 

16:37 

Solar  zenith  angle  (degree) 

34.9 

27.48 

58.6 

44.2 

Cloudiness  {%) 

10.0 

5.0 

30.0 

5.0 

Mean  wind  speed  (m/s) 

6,00 

7.96 

10.0 

5.0 

Average  BP  (1lHg) 

29.92* 

30.26 

29.92* 

29.92* 

Air  mass  type 

1* 

1* 

1* 

1* 

Average  relative  humidity  (%) 

80* 

79.7 

80* 

80* 

Pre  ri  p  liable  water  content  (cm) 

1.5* 

1.5* 

1.5* 

1.5* 

Visibility  (km) 

15.0* 

15.0* 

15.0* 

15.0* 

4 Default  values  used  in  Hydrolight, 


is  ±0,046  rrT1  as  reported  by  the  manufacturer's  specifi¬ 
cations,  Chlorophyll  concentrations  were  determined  from 
ihe  same  samples  using  a  Turner  fluorometer  (Holm- 
Hansen  and  Riemann,  1978).  Details  of  the  instruments 
and  data  processing  of  these  two  stations  are  described  in 
Ivey  et  aL  (2002). 

To  validate  our  approach,  we  need  a  set  of  field  data 
with  biological  measurement  of  Chl-ff(z),  the  above-wa¬ 
ter  measurement  of  /?fJ(A),  and  the  in-water  measurement 
of  lOPs (z).  Apart  from  the  profiles  of  a(zt  A)  and  c(z,  A) 
from  the  ac-9™  absorption  meter  (WET  Labs),  we  espe¬ 
cially  need  the  backscattering  coefficient  bh{z,  A)  for  cal¬ 
culating  the  backseat  ter  fraction  BF(z ,  A).  At  the  time  of 
writing,  there  are  only  two  stations  (LSD52601  and 
LK404U)  with  bb  profiles  measured  from  the  newly  in¬ 
stalled  Hydroscat-6rM  backscattering  meter  (HOB!  Labs). 
Therefore,  we  add  the  data  from  other  two  stations  (Mirirl 
and  Mirir2)  that  have  all  measurements  except  for  bb(zt 
A).  These  measurements  were  obtained  during  a  cruise  in 
the  northern  Gulf  of  Mexico  in  200 L  Details  of  the  in¬ 
struments  and  data  processing  of  these  two  stations  are 
described  in  D'Sa  and  Miller  (2003), 

All  field  data  were  obtained  from  coastal  regions, 
where  LSD52601  and  LK404U  represent  more  clear  wa¬ 
ters  with  moderate  chlorophyll  and  sediment  loading, 
while  Mirirl  and  Mirir2  represent  more  turbid  waters  with 
high  chlorophyll  and  sediment  loading,  A  detailed  descrip¬ 
tion  of  all  stations,  including  their  locations  and  water 
types,  is  given  in  Table  3.  Note  that  the  model  of  Liu  et 
al.  (2002)  can  be  applied  to  Case  I  waters  as  well  as  to 
Case  2  waters  that  is  gelbstoff  rich,  and  the  volume-scat¬ 
tering  phase  function  can  be  generated  dynamically  based 
on  the  backscatter  fraction.  Therefore,  all  field  data  col¬ 


lected  in  these  four  stations  are  valid  inputs  for  the  model 
of  Liu  et  at.  (2002)  to  estimate  the  underwater  light  field. 

3,  Results  and  Discussion 

3.  /  Model-to-data  comparison 

Here  we  present  the  results  of  model-to-data  com¬ 
parison  at  all  stations.  Following  the  Steps  I  and  2  illus¬ 
trated  in  Fig.  2,  the  spectrum  information  of  Rrs(k)  is  sub¬ 
stituted  into  QAA  to  retrieve  Chi ,  Ft  and  BFp .  Since  /?rj(A) 
is  a  weighted  average  signal  over  depths,  the  mean  (hol¬ 
low  symbols)  as  well  as  the  standard  deviation  (error  bars) 
of  the  in  situ  measurements  are  calculated  over  the  e-fold¬ 
ing  depth  and  compared  to  the  retrievals  of  QAA  (filled 
symbols),  respectively  (Fig.  4).  The  e-folding  depth  is 
defined  as  the  depth  where  E0J>AH{z)  has  decreased  to  l/e 
of  its  surface  value.  Note  that  we  compare  BFp  al  band 
488  nm  because  it  is  the  first  waveband  that  both  ac-9 
and  Hydroscat-6  have.  We  also  compare  the  total  absorp¬ 
tion  coefficient  at  440  nm,  which  takes  into  account  the 
CDOM  absorption  through  the  value  of  F  (Eq,  (2)).  For 
the  case  of  the  most  turbid  water  at  station  Mirirl  (Fig, 
4(c)),  the  deviations  between  the  QAA  retrievals  and  the 
means  of  in  situ  measurements  taken  within  the  e-fold¬ 
ing  depth  (only  l  ,8  m)  are  apparent.  The  deviation  can  be 
reduced  if  the  means  are  taken  over  more  depths.  How¬ 
ever,  there  is  no  way  to  determine  to  what  depths  the 
means  should  be  taken  simply  based  on  ihe  existing  meas¬ 
urements,  because  /?„(A)  is  a  weighted  average  signal  over 
depths.  Nevertheless,  Fig,  4  demonstrates  that  QAA  pro¬ 
vides  reasonable  retrievals  of  bulk  lOPs  both  for  verti¬ 
cally  well-mixed  waters  (Fig.  4(a))  and  for  those  cases 
with  vertical  structure  (Figs.  4(b),  (c),  (d)).  These 
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Fig.  4.  Model -to-data  comparison  of  lOPs.  (a)  LSD52601,  (b)  LK404U.  (c)  Mirirl  and  (d)  Mirir2. 


10 


Fig.  5.  Model-to-daia  comparison  of  the  underwater  light  field,  (a)  LSD5260L  (b)  LK404U,  (c)  Mirirl  and  (d)  Mirir2* 


retrievals  are  used  as  inputs  to  the  model  of  Liu  et  aL 
(2002)  to  calculate  E0JfAk(z),  which  corresponds  to  the  in 
situ  measurements  well  (Fig.  5), 

As  mentioned,  light  absorption  by  phytoplankton  or 
other  water  constituents  generally  have  no  preferential 
direction  (Kirk,  1 994)*  Therefore,  we  should  take  the  sca¬ 
lar  irradiance  E0  rather  than  the  downwelting  planar  irra- 
dtance  Ed  into  account.  Because  the  in  situ  measurements 
are  Ejfar{i)  which  cannot  be  compared  directly  to  our 
model  output  EQPAR(z),  we  use  H42  as  a  surrogate  to  simu¬ 
late  both  E0tFAR{z)  and  EdPAR{z)  under  the  same  computa¬ 
tional  conditions.  Figure  5  shows  that,  using  the  retrievals 
of  Chi ,  F  and  BFp  from  QA  A,  H42  offers  a  simulation  of 
EdPAR(z)  that  corresponds  closely  to  the  in  situ  measure¬ 
ments.  In  the  meantime,  the  results  from  the  model  of 
Liu  et  aL  (2002)  are  virtually  identical  to  H42  simulations 
of  Eq  par(z).  This  comparison  verifies  our  approach. 

A  major  premise  of  this  work  is  the  vertical  homo¬ 
geneity  of  IOPs,  which  is  not  always  valid  (e.g.  Figs,  4(b) 
and  (c))*  There  is  currently  no  practical  way  to  resolve 
the  vertical  structure  simply  from  surface  observations 


of  Rrs.  The  spectrum  of  Rr ,  should  be  signals  weighted  by 
KJ^X)  over  the  first  optical  attenuation  length,  dependent 
on  wavelength  and  the  optical  properties  of  the  water  at 
that  wavelength  (Platt  and  Sathyendranath,  1988),  Nev¬ 
ertheless,  both  the  GC2-v4  algorithm  and  our  approach 
provide  good  retrievals  that  represent  the  bulk  properties 
in  the  upper  ocean.  The  model-to-data  comparison  dem¬ 
onstrates  that  a  good  estimate  of  EQ  PAR(z)  at  depth  can  be 
obtained,  as  long  as  the  bulk  values  of  IOP  are  retrieved 
for  the  important  surface  layer*  A  detailed  distribution  of 
the  IOPs  in  the  vertical  direction  would  be  of  great  help 
in  estimating  EQPAR(z)  at  depths*  However,  such  informa¬ 
tion  usually  derives  from  simulations  of  biogeoc  hem  teal 
models  rather  than  remote  sensing  (e.g.  see  Bissett  et  aLt 
1999).  The  model  of  Liu  et  aL  (2002)  needs  to  be  further 
expanded  to  deal  with  a  stratified  water  column,  if  the 
goal  is  to  incorporate  the  optical  model  into  a 
biogeochemical  model. 

Figure  5  also  illustrates  the  difficulties  in  obtaining 
a  smooth  profile  of  EdPAR{z)  from  in  situ  measurements. 
The  irradiance  fluctuation  is  a  typical  phenomenon  that 
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Fig*  6,  Model -to -mode  l  comparison  of  E$pAR{z)l  where  X  rep¬ 
resents  various  approaches,  including  H42  { Hydro! ight  ver¬ 
sion  4.2),  K£j( 490)  (Mueller,  2000),  KEiI(A)  (Eq.  (6)),  and 
Liu  et  al.  (this  work). 


is  not  due  to  small-scale  vertical  variations  in  the  IOPs 
but  almost  entirely  due  to  the  wavy  surface  (Zaneveld  et 
aL%  2001),  Another  limitation  of  the  model-to-data  com¬ 
parison  is  that  the  dataset  collected  from  in  situ  measure¬ 
ments  only  covers  a  limited  range  of  IOPs.  To  further 
verify  our  approach  for  a  large  variety  of  water  types,  we 
conducted  the  following  model-to-model  comparisons. 

3.2  Model  -  to-  m  ode  I  comparison 

The  procedure  used  for  model -to-model  comparisons 
is  briefly  as  follows.  For  a  set  of  computational  condi¬ 
tions  of  IOPs  and  the  ambient  optical  environment,  the 
H42  radiative  transfer  model  (Mobley  and  Sundman, 
2001)  is  used  to  simulate  £H,(490),  £h,{555)  and  Rrs( A).  A 
detailed  vertical  profile  of  E^pAR{z)  is  also  generated  to 
represent  the  simulated  underwater  light  field  derived 
from  a  full  solution  of  the  radiative  transfer  equation.  The 
latest  formula  (Mueller,  2000)  is  used  to  calculate 
***490)  from  £w(490)  and  Lw( 555),  The  ratio  of  EJEj  is 
nominally  l, 4- 1.8  up  to  2. 0^2,5  in  very  turbid  waters 
(Kirk,  1994).  Therefore,  A^490)  is  assumed  to  approxi¬ 
mate  K0  PAH  averaged  over  the  entire  PAR  spectral  range 
and  depth,  and  (z)  is  estimated  by  attenuating 

the  below-surface  value  of  EQPAR( 0‘)  by  KEJ( 490), 
Austin  and  Petzold  (1986)  found  a  relationship: 


^ Ed  W  ~  wafer  W  ]Lf(X\ 

^  Ed  (490)  -  KEd  wntcr  ( 490)  1 


(6) 


which  allows  us  to  estimate  KEJ,K)  from  ^^(490).  There¬ 


Fig.  7.  Model-to-modd  comparison  of  the  relative  errors  e(z) 
in  predicting  E0  Par{z)  by  use  of  (a)  KEt^X)  approach,  and 
(b)  Liu  et  al.  approach.  The  controlled  variable  $s  is  varied 
from  0*  to  60°. 


fore,  another  profile  E^pj^(z)  can  be  obtained  by  at¬ 
tenuating  the  below-surface  value  of  EQJ>AR( 0“)  by  fCEd(X). 
Finally,  the  H42  simulated  /?fJ(A)  and  EjlK  0*)  curves 
are  used  to  calculate  E^pARlt{z)  following  the  approach 
illustrated  in  Fig,  2, 

We  start  from  a  simulation  of  a  typical  Case  l  water 
and  compare  the  predictions  of  E0J>AR(z)  using  different 
approaches  (Fig.  6).  A  detailed  description  of  this  set  of 
computational  conditions  is  given  in  the  caption  of  Fig. 

6,  The  relative  errors  of  are  as  high  as  188.0% 

at  the  euphotic  depth  (38.0  m)t  while  £^j*a/(z)  closely 
matches  the  full  H42  solution  E^ARAR{z)  (8.1  %  at  the  same 
depth),  EQ^R  f  (z)  indeed  gives  a  better  prediction  than 

(2)  by  taking  the  spectral  variation  of  K£j[X)  into 
consideration.  However,  apparent  deviations  between 
Eq,pa^(z)  an^  can  still  be  seen  from  5  m  to  30 
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Fig.  8.  As  Fig.  7  except  that  the  controlled  variable  F  is  varied 
from  0.1  to  1,6, 


m.  The  relative  error  of  E*p}^(z)  reaches  a  maximum 
value  of  17,0%  at  the  depth  of  14.0  m,  while  the  relative 
error  of  is  merely  2.1%  at  the  same  depth. 

This  error,  as  mentioned  above,  resulted  from  neglecting 
the  variation  of  KEJ  with  depth.  To  validate  our  approach 
and  assess  the  error  that  might  result  from  other  sources, 
we  extend  the  model-to-model  comparison  by  varying  one 
computational  condition  at  a  time  and  keeping  the  others 
the  same  as  those  used  in  Fig,  6. 

Figure  7  shows  the  result  of  varying  the  solar  zenith 
angle  8S  from  0°  to  60°,  which  can  be  used  to  illustrate 
the  variation  of  the  underwater  light  field  within  one  day. 
Even  for  a  typical  Case  I  water,  0S  has  a  significant  ef¬ 
fect  on  E0PAR(z)  that  cannot  be  modeled  accurately  using 
the  Kf/h)  approach.  The  relative  errors  range  from 
-30,0%  to  54.9%  at  different  depths  (Fig.  7(a)),  By  con¬ 
trast,  our  approach  gives  an  estimation  that  reduces  the 
relative  error  (0  to  16.3%),  regardless  of  the  different 
used  (Fig,  7(b)).  Varying  the  solar  zenith  angle  $s  would 
have  similar  effects  on  the  simulation  of  underwater  light 
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Fig.  9.  As  Fig.  7  except  that  the  controlled  variable  BF  is  var¬ 
ied  from  0.0123  to  0,0363. 


field  to  varying  other  parameters  of  incident  light  field, 
such  as  the  surface  wind  speed,  visibility,  cloudiness,  etc. 

Figure  8  gives  the  result  of  changing  the  CDOM  ra¬ 
tio  F  from  0.1  to  L6,  The  relative  errors  range  from 
-58.4%  to  5,2%  at  different  depths  (Fig.  8(a)).  By  con¬ 
trast,  our  approach  gives  an  estimation  that  reduces  the 
relative  error  (-10.3%  to  13.9%),  regardless  of  the  dif¬ 
ferent  F  used  (Fig.  8(b)).  Figure  9  shows  the  result  of 
changing  the  backscatter  fraction  BFp  from  0.0123  to 
0.0363,  The  relative  errors  range  from  -19.0%  to  66,9% 
at  different  depths  (Fig.  9(a)).  By  contrast,  our  approach 
gives  an  estimation  that  reduces  the  relative  error  (-0.8% 
to  1 1 ,4%),  regardless  of  the  different  BFp  used  (Fig.  9(b)), 
The  KeJ^X)  approach  (Eqs.  (1)  and  (6))  is  developed  for 
Case  1  waters.  If  the  lOPs  are  deviated  from  the  assump¬ 
tion  of  Case  1  waters,  it  is  expected  that  using  the  Kf/X) 
approach  to  predict  E0  f,AR(z)  will  generate  a  significant 
error  (Figs.  8(a)  and  9(a)).  By  contrast*  our  approach  pro¬ 
vides  an  estimation  that  closely  matches  the  full  H42  so¬ 
lution,  regardless  of  the  different  values  of  F  or  BFp  used 
(Figs.  8(b)  and  9(b)).  Varying  the  CDOM  ratio  F  or  the 
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Fig.  10.  Applications  in  processing  ScaWiFS  imagery  taken  during  the  EcoHAB  cruises  in  2001  on  the  West  Florida  Shelf, 
(a)  Chlorophyll -a  concentration  (SeaWiFS  OC2V4  algorithm),  (b)  Chlorophyll-a  concentration  (Liu  et  at.),  (c)  a(443)  (Liu  ei 
alX  (d)  />^(443)  (Liu  et  al>),  (e)  Chlorophyll-^  concentration  (Liu  et  al and  (f)  Dc  (Liu  et  o/t). 
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backscatter  fraction  BFp  would  have  similar  effects  on 
the  simulation  of  underwater  light  field  to  varying  other 
parameters  of  water  constituents,  such  as  the  chlorophyll 
concentration, 

4.  Processing  SeaWiFS  Imagery 

The  newly-developed  approach  was  applied  to  proc¬ 
ess  SeaWiFS  imagery  taken  during  the  EcoHAB  cruises 
in  2001  on  the  West  Florida  Shelf.  Figure  10(a)  shows 
the  distribution  of  chlorophyll-^  concentration  derived 
from  the  standard  SeaWiFS  OC2V4  algorithm.  The  white 
areas  indicate  those  pixels  that  are  either  masked  by  cloud 
or  failed  the  quality  lest,  while  the  white  cross  mark  pin¬ 
points  the  place  where  the  in  situ  measurements  (Figs, 
4(b)  and  5(b))  were  made.  Since  our  approach  only  re¬ 
quires  o4[2,  a443  and  hbp443  (see  Step  2  in  Fig.  2),  the 
retrievals  of  these  values  by  use  of  QAA  are  shown  in 
Figs,  10(b),  (c)  and  (d),  respectively.  As  expected,  the 
distribution  illustrates  the  influences  of  coastal  runoff  on 
the  IOPs.  The  retrieval  of  chlorophyll-a  concentration 
following  Step  2  of  Fig,  2  is  illustrated  in  Fig.  10(e). 
Compared  to  the  results  of  OC2V4,  our  approach  gives 
comparable  values  of  Chi  in  general.  Figure  10(0  shows 
the  depth  of  the  euphotic  zone  where  EQ  PAR  equals  1%  of 
the  subsurface  value  of  EQ^AR.  This  depth  is  usually  treated 
as  an  approximation  of  the  compensation  depth  Dc  (Mann 
and  Lazier,  1996),  which  plays  a  significant  role  in  the 
ecosystem  of  phytoplankton.  Failure  to  correctly  estimate 
the  underwater  light  field,  and  hence  the  compensation 
depth,  will  generate  significant  errors  in  modeling  the 
plankton  ecosystem, 

5,  Conclusions 

Regional  to  global  estimates  of  the  underwater  light 
field  are  required  to  develop  realistic  models  of  ocean 
primary  productivity,  and  heat  and  carbon  fluxes.  To  date, 
most  attempts  to  incorporate  satellite  ocean-color  data  are 
based  on  a  simple,  but  error-prone  ^^490)  model.  We 
present  a  new  approach  that  incorporates  two  models  (Lee 
ei  a/.,  2002;  Liu  et  ai. ,  2002)  to  estimate  the  underwater 
light  field  from  remote  sensing  of  ocean  color.  A  thor¬ 
ough  model-io-data  and  model-to-model  comparison  was 
performed  to  verify  this  new  approach.  We  successfully 
mapped  the  compensation  depth  by  applying  this  new 
approach  to  process  SeaWiFS  imagery.  The  current  model 
of  EQ  PAJi(z)  is  limited  to  a  vertically  homogeneous  water 
column  with  variables  in  certain  ranges:  A  (400-700  nm), 
Chi  (0-10  mg  nT3)  and  BFp  (0.01-0.04)  (Liu  etaL ,  2002). 
Future  work  will  improve  the  model  by  removing  some 
of  the  limitations,  and  develop  a  better  scheme  of  inver¬ 
sion  to  retrieve  IOPs.  This  work  indicates  that  incorpo¬ 
rating  two  models  (Lee  et  aLy  2002;  Liu  et  ai.t  2002)  may 
provide  large-scale  estimates  of  the  underwater  light  field 
that,  while  not  perfect,  are  significantly  more  accurate 


than  a  Ked  approach  and  computationally  feasible  for  the 
generation  of  routine  products.  The  product  of  EQFAM(2) 
can  be  obtained  routinely  from  ocean-color  data  and  may 
have  significant  implications  for  the  estimation  of  global 
heat  and  carbon  budget. 
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Appendices 

Symbol 

Units 

Description 

a 

m_l 

total  absorption  coefficient 

a;  (A) 

normalized  chlorophyll-specific  ab¬ 
sorption  coefficient 

m_l 

absorption  coefficient  of  CDOM 

nr1 

absorption  coefficient  of  pure  water 

parameter  required  for  McCormick 
model  of  JJ  (2) 

parameter  required  for  McCormick 
model  of  Ji  (z) 

*2 

parameter  required  for  McCormick 
model  of  /j  (2) 

BF 

backscatter  fraction 

BFo 

backscatter  fraction  of  particles 

BFW 

backscatter  fraction  of  water 

b 

HI-1 

total  scattering  coefficient 

h 

nT1 

backscattering  coefficient 

bp 

m"1 

scattering  coefficient  of  particles 

K 

nT1 

scattering  coefficient  of  pure  water 

Chi 

mg  nT3 

chlorophyll-a  concentration 

Oc 

m 

compensation  depth 

Bo 

WmV1 

scalar  irradiancc 

Eft. PAR 

Wm‘2sr‘l 

wavelength-integrated  scalar  irradi- 
ance  in  the  PAR  range 

Ed 

Wm"2sr“f 

downwelling  planar  irradiancc 

Bjpar 

Wm'2sr_l 

wavelength-integrated  downwelling 
planar  irradiancc  in  the  PAR  range 

F 

CDOM  ratio 

L 

Wm^sr^nm" 

1  radiance 

Lw 

Wm'^r’nm' 

1  water-leaving  radiance 

Eej 

nT1 

total  diffuse  attenuation  coefficient  for 
downwelling  planar  irTadiance 

^£t/.wofEr 

m"1 

same  as  KEd  but  for  pure  water 

F 

parameter  required  for  McCormick 
model  of  Ji(z) 

Q 

parameter  required  for  McCormick 
model  of  ji  (2) 

Ks 

the  remote  sensing  reflectance  mea¬ 
sured  just  above  the  surface 

rn 

the  remote  sensing  reflectance  mea¬ 
sured  just  below  the  surface 
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normalized  phase  scattering  function 
normalized  phase  scattering  function  of 
pure  water 

normalized  phase  scattering  function  of 
particles 
wavelength 
average  cosine 
solar  zenith  angle 
scattering  angle 
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